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The random sequential adsorption of various particle shapes is studied in order to determine 
the influence of particle anisotropy on the saturated random packing. For all tested particles 
there is an optimal level of anisotropy which maximizes the saturated packing fraction. It is found 
that a concave shape derived from a dimer of disks gives a packing fraction of 0.5833, which is 
comparable to the maximum packing fraction of ellipsoids and spherocylinders and higher than 
other studied shapes. Discussion why this shape is beneficial for random sequential adsorption 
is given. 


1 Introduction 

Random sequential adsorption (RSA) is a fundamental problem 
that has attracted unabated interest for decades. In 1939, Flory 1 
studied the attachment of blocking pendant groups on a lin¬ 
ear polymer, which is effectively a one-dimensional RSA prob¬ 
lem. Renyi 2 introduced another famous one-dimensional RSA 
problem—the parking of cars along an unmarked curb. Feder 3 
helped to make RSA a very popular tool for modeling monolayers 
obtained as a result of irreversible adsorption 4 ®. More recently, 
random packings generated by RSA have been of interest in a 
number of scientific fields, e.g., soft matter 7 9 , surface science®], 
mathematics^ 11 , telecommunication 12 and information theory 13 . 

The RSA algorithm is based on consecutive tries to add a par¬ 
ticle to a packing. Firstly, a particle’s position (and orientation in 
case of anisotropic objects) is drawn according to the probability 
distribution which reflects the structure of an underlying surface 
- a homogeneous surface corresponds to the uniform probability 
distribution. Then the particle is tested if it overlaps or intersects 
with any of particles already added to the packing. If not, the 
particle is added to the packing; otherwise it is abandoned. These 
steps should be continued until the packing is saturated, i.e., there 
is no room for any additional particle on a surface. However, typ¬ 
ically an algorithm is stopped when the probability of successful 
adding of a particle is sufficiently small, and extrapolations are 
made to estimate the maximum “jamming" coverage. 

The properties of packings generated by RSA algorithm 


have been checked for a number of different particle shapes, 
e.g., spheres ® 14 ! 15 !, spherocylinders and ellipsoids 1617 , rectan- 
glesESHS] anc [ polymers®! 21 . Results obtained for anisotropic 
particles show that saturated random coverage fraction reaches 
its maximum for moderate anisotropy, i.e., when long-to-short 
particle axis ratio is approximately 1.5-2.0® 19 . The highest sat¬ 
urated coverage fraction obtained was 0 max = 0.583 ±0.001, for 
both ellipsoids and spherocylindersGSE 17 . Both of these are con¬ 
vex shapes. The highest saturated coverage fraction for concave 
particles was observed for a dimer built of two overlapping disks, 
with 0 max = 0.5793 ±0.0001®!, slightly lower than above values. 

A similar problem of finding the highest possible random close 
packing of various objects in three dimensions is crucial for under¬ 
standing granular and soft matter properties.®!. However close 
packings differs significantly from saturated packings generated 
by the RSA algorithm as in the first case particles are tightly 
packed and are touching most of their closest neighbors. It is 
interesting that the random close packing fraction of ellipsoids 
can be as much as 0.74 while for spheres it is 0.64®®. More¬ 
over, it occurs that the highest packing were found for particles 
(ellipsoids, spherocylinders or dimers) of long-to-short axis ra¬ 
tio around 1.5, which is similar as for RSA in two dimensions. 
These results were confirmed both: numerically®®! 27 and theo¬ 
retically®]. 
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The primary aim of this study is to analyze set of different con¬ 
cave shapes in order to find if one of them can give a higher sat¬ 
urated random packing fraction than previously found. This will 
help in designing particles, such as nanoparticles, that can cover a 
surface most efficiently. Additionally, we study the kinetics of RSA 
for nearly symmetric particles as they appear to be very sensitive 
on particle shape anisotropy®®!. 
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2 Model 

The shapes we consider are shown in Fig. [l] Dimers, trimers, 


«tit# 

Fig. 1 The five types of shapes for which saturated random packings 
are studied. Disks have a unit radius, and parameter x e [0,1] 
corresponds to half the distance between closest spheres or squares. 
Squares have a unit side size, thus parameter a equals y/2/2. 


tetramers and hexamers are built of identical disks of unit radius, 
and the distance between neighboring disks centers is 2x. The 
parameter v if not stated otherwise was taken from the interval 
[0,1]. For jc = 0 all these shapes are equal to a single disk, which 
of course is isotropic, while for x = 1 each disk touches two neigh¬ 
boring disks as presented in the upper row of Fig. [l] Note that the 
space inside ring of disks also belongs to the particle. Such choice 
of shapes originates fronJSl, where properties of dimer random 
packings were carefully studied. Here we want to test how the 
more complex shape of particle shape affects properties of sat¬ 
urated random packing. The last shape is built of two partially 
overlapping squares of a unit side size. Parameter a shown in 
Fig. [l] equals to y/2/2. This shape is an analogue of generalized 
dimer but instead of disks, it is built of two squares. The sur¬ 
face area, for 0 < x < 1 , of dimer GS 2 ), trimer (S 3 ), tetramer (S 4 ), 
hexamer (£ 5 ) and the square dimer GS^s) presented in Fig. |T] are 
given respectively: 

S 2 (x) — 2tt — 2 ^arccosx — xy/Y^-x 2 ^ (1) 


+ 2 arcsin ^ ^ \/ 3 — 3 x 2 ^j 

± 2 arcsin ^ - (v + y/3 — 3 x 2 ^j 

± ijc\/ 2 jc 2 — 2jc \/3 — 3jc 2 + 1 

+ \^\J~2j{ 2 4-2 

± ^- \j [x 2 ~ l) (—2x 2 + 2x\f3~— 3x 2 — ij 
— yj(x 2 — l) 2x 2 — 2xy /3 — lx 2 — 1 ^ (4) 

S 2 s(x) = 2(\+2x-x 2 y, (5) 

where in S 3 (x) and S^(x), the upper sign in ± or =p is for 0 < v < 
1 / 2 , and the lower sign is for 1 /2 < v < 1 . 

These shapes were thrown onto a square surface of a side size 
up to 1000 with an area up to Sc = 10 6 . We decided to use open 
boundary conditions similar tJ2l. Such a choice allowed us to 
use data published there for comparison purposes. The packing 
fraction is equal to: 

e(t)=N(t)^= P (t)S p , ( 6 ) 

Sc 

where S p is a particle surface area and N(t) is a number of parti¬ 
cles in a packing after a number of RSA iterations corresponding 
to time t measured in the dimensionless time units 



where n is number of RSA algorithm steps, and p(t) = N(t)/Sc is 
the surface density of particles in a packing. 

The simulation was stopped when t = 10 5 . At this moment, the 
probability that a randomly chosen place is large enough for an 
additional particle is well below 10 -6 . The number of particles 
in a single random packing was at the order of 10 5 . To improve 
statistics, up to 100 independent simulations were performed for 
each shape. During the simulations the number of particles in a 
packing N(t) was recorded. 


S 3 (x) = tt(3 =f 1 )/2 + arcsin x + y/3x 2 -\-4xy/ 1—x 2 


+ 

arcsin ^ (^c — y/3 — 3 x 2 ^j ^ 


± 

arcsin ^ + y/3 — 3x 2 ^j ^ 


± 

^x\l 2x 2 — 2 xyj 3 — 3x 2 + 1 


+ 

y>c\l 2x 2 + 2xy/3 — 3x 2 + 1 


± 

^ (x 2 — l) (^—2x 2 + 2xy/3 — 3x 2 — 1 j 


- 

yj(x 2 — l) 2x 2 — 2xy/3 — 3x 2 — 1^ 

( 2 ) 

s 4 0 ) 

= 7T + 4x 2 + 4xy/1 — x 2 + 4 arcsinv 

(3) 


3 Results 

Fragments of sample packings are presented in Fig. [ 2 ] Examples 
of dimer packings have been shown in 22 . The total number of 
RSA algorithm steps is large, but finite. Thus, the resulting pack¬ 
ings may be not be saturated. To find number of particles in sat¬ 
urated packing from a finite-time simulation a knowledge about 
RSA kinetics is essential. Therefore, to study saturated packing 
fractions the kinetics of the process has to be clarified first. 

3.1 RSA kinetics 

The RSA kinetics for large enough time t is governed by the power 
law^ 

0(t) = 0 max -Ar l / d ( 8 ) 


S&(x) = tt(3 =f 3)/2 + 2 arcsin x + 6y/3x 2 + %xy /1 — x : 

+ (1 ± 1) j^arcsin ^y/3xj + arcsin ^y/ 1 — 3x 2 ^j j 


where 0 max = d(t °°), A is a positive constant and d depends 
on particle shape and properties of a surface. For flat and homo¬ 
geneous surfaces d can be interpreted as a number of degrees of 
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ior. In other words, the apparent change of d from 3 to 2 as i 
decreases represents crossover behavior between the two expo¬ 
nents. For trimers this drop occurs for very small x, which resem¬ 
bles its behavior for dimersEH. In the case of square dimers very 
slight decrease is observed, because even for x = 0 the particle 
is still anisotropic. This is in agreement with previous observa- 








Fig. 2 Examples of obtained random packings for studied shapes after 
t = 10 5 iterations of RSA algorithm. Parameter x equals 0.5 for all four 
shapes. At this point the systems were at approximately 99% of the 
extrapolated maximum packing. 


freedom of a particleESE. Thus, for RSA of disks on two dimen¬ 
sional flat surface d = 2 but for anisotropic particles d = 3 as an 
orientation of a particle is an additional degree of freedom. It is 
worth noting, that even when particle anisotropy is very small, 
d = 3 describes the ultimate asymptotic behavioJ 16 l 22 l 29 [ 

The RSA kinetics of our shapes are presented in Fig. [3] Data 



3.2 Saturated random packing fractions 

Knowing parameter d and substituting y = t~ l / d , Eq. ^ can be 
rewritten in a form: N(y) = A^max+A'y. Thus, points (N(y),y) mea¬ 
sured during an RSA simulation lie along a straight line which 
crosses the axis y = 0 at A max . It solves the problem of finite time 
simulations. The error of such A max estimation originates in er¬ 
ror of the exponent — 1/d, which in our simulation is below 0.005 
(see Fig|3]). The corresponding component to the relative error of 
Nmax is 0.06% (see Fig|4|). Its value is comparable with the statisti- 



Fig. 4 (Color online) Estimation of number of particles in a saturated 
random packing. The main plot shows influence of 1/d exponent error 
on obtained number of particles in saturated packing N max for dimer 
characterized by x = 0.5. Dots represent data taken from simulations for 
y = t~ l / d±0 005 . Lines are linear fits: N(y) = 1.1464-10 5 -48598y and 
N(y) = 1.1456-10 5 —51573y. Inset shows dependence of N max on surface 
size. Dots represent data taken from simulations while line is a quadratic 
fit: N = 162.64-0.70174^ + 0.114995. 




Fig. 3 (Color online) RSA kinetics for trimers, tetramers, hexamers and 
square dimers. Main panels show number of added particles in a time 
unit versus time. Insets show fitted exponent in Eq. [^dependence on x. 
Error bars are smaller than symbol sizes. Dashed lines correspond to 
the d = 3 characteristic behavior for anisotropic molecules. 

for dimers has been presented inEU. For all shapes Eq. ^ is 
fulfilled as the data obtained from numerical simulations lays 
along straight lines in a log-log scale. The parameter d obtained 
from fitting numerical data to the relation © is, as expected for 
anisotropic particles, around 3. In a limit of small x, for tetramers 
and hexamers the drop of d down to 2 is observed, which is ex¬ 
pected, as the particle shape arises from disk-like geometry and 
the simulations have not yet reached the true asymptotic behav- 


cal error coming from averaging over 100 different simulations. It 
is worth mentioning that the error component coming from finite 
time simulations can be totally eliminated for spherical particles 
using an algorithm presented inR 

which allows one to obtain 
saturated packings in a finite-time simulations by tracing regions 
where addition of particles is possible. Thus, the RSA sampling 
can be limited to these regions only. Unfortunately, in the case of 
anisotropic particles, the complexity of this algorithm raises as the 
success of placing a particle depends not only of its position but 
also on its orientation. Another problem originates in the finite 
size of a system, which is especially important as periodic bound¬ 
ary conditions are not used here. It is expected that in general the 
number of particles in a packing are a quadratic function of the 
surface size: N{Sq) = aSc + byfSc + c (see the inset in Fig|4]). Thus 
for an infinite system (, Sc -+ °°), a packing density N(Sc)/Sc = a. 
Therefore, to get a packing density in a limit infinitely large sys¬ 
tems, several different sized systems should be analyzed to get 
the a coefficient in the above relation. In this study Sc varies 
from 2.5 • 10 5 up to 10 6 . Results presented in Fig(4] suggest that 
a systematic relative error coming from a finite surface size is in 
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our case at the order of 0.1% and is the biggest among all error 
components. Therefore, the best way to improve the quality of 
obtained results is to simulate RSA on larger surfaces. 

Obtained packing fractions are presented in Fig. [5] For all 




X 


Fig. 5 (Color online) Saturated packing fraction dependence on 
parameter x for trimers, tetramers, hexamers and square dimers. Dots 
represents data from numerical simulations. Error bars are smaller than 
symbol sizes. Solid red lines are quadratic fits near maxima of 0 max . 
Data for dimers were taken from^j. The fits are: 

0 2 = 0.5483 + 0.12089x - 0.1175 lx 2 , 

0 3 =0.54933+ 0.088299x-0.091899x 2 , 

0 4 = 0.54835+ 0.033967X-0.050194X 2 , 

0 6 = 0.54867+ 0.017815x —0.042357X 2 , 

0 2 , =0.52839 + 0.10979x-0.37802x 2 . 


shapes the maximal packing fraction is reached for x e [0,1]. This 
confirms reasoning presented in® ' that anisotropy prefers par¬ 
allel alignment of particles for large t, which causes increase of 
packing fraction, but, on the other hand, at the beginning of RSA 
the anisotropic particle blocks significantly more space than its 
surface area, which lowers the packing fraction. Thus, the opti¬ 
mum is reached for a small anisotropy. Fitting a quadratic func¬ 
tion to the data near the maximum allows us to estimate an opti¬ 
mal anisotropy, which in our case is determined by parameter x, 
and the value of the highest possible packing fraction. For conve¬ 
nience these data are collected together in the Table [l] 

4 Discussion 


Table 1 Maximal possible saturated packing fractions and 
corresponding values of parameter x for which they are reached. The 
statistical error of 0 max does not exceed 0.0005. 


shape 

X 

^max 

dimeiH 

0.5098 

0.5793 

trimer 

0.4804 

0.5705 

tetramer 

0.3384 

0.5541 

hexamer 

0.2103 

0.5505 

square dimer 

0.1452 

0.5364 


spherocylinders (0.583) reported in 16 . However, these dimers, 
trimers tetramers and hexamers do not have the best possible 
shape to reach highest random packing fraction. As illustrated in 
Fig. [6] the small blue areas in both sides of a dimer cannot be oc- 




Fig. 6 (Color online) Top panel: the transition from generalized dimer to 
a shape with larger surface area and the same blocking area. Bottom 
panels: packings at the end of simulation bulid of such shapes for 
x = 0.6 (left) and x = 1.2 (right). 


cupied by any other particle. It means that this surface is blocked 
but it is not counted in the packing fraction. There is a direct 
projection between packing composed of dimers and rightmost 
particles in Fig. [6] Therefore saturated random packing fraction 
for this shape will be slightly higher than for dimers. Similar rea¬ 
soning can be performed for trimers, tetramers and hexamers. 
Note that this reasoning is valid also for x > 1, however, too large 
a value of x can destroy the mentioned projection between two 
types of shapes. For example for dimer the projection is valid for 

X < y/2. 

The surface area of one of these additional regions is 


^ad M 


X ^a/ 4 — x 2 — y/\ — X 2 ^ 
x V4 — x 2 - § 


— arcsinx 


0 < x < 1 
1 <x< y/3 


(9) 


For all shapes we studied, the maximal possible packing fraction is 
smaller than for generalized dimers®. The more complex shape 
type, the lower the value of the maximum packing fraction. The 
lowest values are obtain for square dimer. In this last case the 
0 m ax is smaller than for rectangles (0. 55)^5® as well as for disks 
(0.5470) 15 . In all cases it is smaller than value for ellipsoids or 


Taking this into account and performing the same analysis (see 
Fig. 0 as in Sec. |3.2| the following values are obtained (see Ta¬ 
ble |2|. Note that in case of dimers with x > 1 the additional sur¬ 
face does not compensate a smaller number of particles in the 
packing so the packing fraction lowers with growing x. There¬ 
fore we focused on x < 1 in the case of trimers, tetramers and 
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Fig. 7 (Color online) Saturated packing fraction dependence on 
parameter x for dimers, trimers, tetramers and hexamers with additional 
blocked regions taken into account (“Smoothed" shapes). Dots 
represents data from numerical simulations. Error bars are smaller than 
symbol sizes. Solid red lines are quadratic fits near maxima of 0 max . 
Data for dimers were taken frorr® The fits are 
0 2 = 0.55513 + 0.089974X-0.072349X 2 , 

0 3 = 0.55253+ 0.070875x-0.057558X 2 , 

0 4 = 0.55062 + 0.021978x-0.022673x 2 , 

0 6 =0.54883+ 0.015239X-0.024414X 2 . 


hexamers. The calculated packing fractions are bigger, but the 

Table 2 Maximal possible saturated packing fractions and 
corresponding values of parameter x for which they are reached. Here 
the additional blocked regions were taken into account (“smoothed” 
shapes). The statistical error of 0 max does not exceed 0.0005. 


shape 

X 

$max 

dimer 

0.6347 

0.5833 

trimer 

0.6146 

0.5744 

tetramer 

0.4969 

0.5560 

hexamer 

0.3036 

0.5513 


difference is less than 0.005 and becomes smaller with growing 
number of disks in particle. On the other hand, the value ob¬ 
tained for dimeric particle is now similar to the one for ellipsoids 
and spherocylinders. As expected, the value of x where the pack¬ 
ing reaches a maximum is larger when we include the blocked 
areas in the particles. For the case of the smoothed out dimer, we 
are finding a maximum coverage at x = 0.6347, or in other words 
the aspect ratio a (ratio of maximum dimensions in the two direc¬ 
tions) is 1.6347. This compares with a maximum for other shapes 
as given in Table [ 3 ] Evidently, by making the smoothed concave 
dimer shape, one is able to reach a coverage comparable to the 
maximum of certain convex shapes (the ellipse and spherocylin- 
der). This is because in RSA there is a tendency for particles of 
moderate aspect ratio to adsorb efficiently in both parallel and 
perpendicular directions. Having an indentation in the sides of 
the particles allows one of the disk ends of a neighboring particle 
to occasionally fit in nicely, and thus increasing the overall cov¬ 
erage. With longer aspect ratios, large open spaces form which 
lowers to overall coverage. 

The quadratic fits given in captions of Figs. [ 3 ] and [ 7 ] are valid 


Table 3 Comparison of the maximal coverages for various oblong 
shapes, a is the aspect ratio, equalling 1 +x for dimers. Numbers in 
parentheses are errors in the last digit. 


shape 


rectangle 

dimer 

ellipse 

spherocylinder 
smoothed dimer 


a 

$max 

Ref. 

1.618 

0.553(1) 

Eg] 

1.5098 

0.5793(1) 

m 

2.0 

0.583(1) 

G! 

1.75 

0.583(1) 

Hi 

1.6347 

0.5833(5) 

this work 


only in a close neighborhood of the maximum coverage. If some¬ 
one is interested in finding saturated packing fraction of dimers 
for arbitrary x it is better to use the p (x) dependence presented in 
the caption of Fig. [ 3 ] The saturated packing fraction is then given 



Fig. 8 (Color online) The dependence of the surface density of dimers 
(both standard and smoothed) on parameter x. Dots represents data 
from numerical simulations. Error bars are smaller than symbol sizes. 
Solid red lines is a 4 th order fit of 

p(x)=y = 0.17452-0.17742x + 0.14968x 2 -0.079263x 3 + 0.019315x 4 . 


by: 


0max(X) — S( x )P( x )i (10) 

where S(x) = S 2 (x) for standard dimer and S(x) = S 2 (x) -\-2S^(x) 
for smoothed one. 


5 Summary 

Several anisotropic and concave shapes of particles were ana¬ 
lyzed in terms of maximal possible random packing fraction. It 
was found that the highest packing fraction is obtained for parti¬ 
cle built of two disks with additional blocking area (see Fig. [6]), 
where the distance between their centers is 2 x 0.6347 which in¬ 
terestingly is more than 20% larger than in the case without the 
extended blocking region. The saturated random packing fraction 
for such a shape is 0.5833 ±0.0001 and is at the level of the high¬ 
est of previously reported values for other particle shapes. Future 
studies of higher precision can show which of the high coverage 
shapes (smoothed dimer, ellipse, or spherocylinder or other ob¬ 
long object) is the absolute best. 
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